Heavy sterile neutrinos, entropy and relativistic energy production, and the relic 

neutrino background 
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We explore the implications of the existence of heavy neutral fermions {e.g., sterile neutrinos) 
for the thermal history of the early universe. In particular, we consider sterile neutrinos with rest 
masses in the 100 MeV to 500 MeV range, with couplings to ordinary active neutrinos large enough 
to guarantee thermal and chemical equilibrium at epochs in the early universe with temperatures 
T > 1 GeV, but in a range to give decay lifetimes from seconds to minutes. Such neutrinos would 
decouple early, with relic densities comparable to those of photons, but decay out of equilibrium, 
with consequent prodigious entropy generation prior to, or during, Big Bang Nucleosynthesis (BBN). 
Most of the ranges of sterile neutrino rest mass and lifetime considered are at odds with Cosmic 
Microwave Background (CMB) limits on the relativistic particle contribution to energy density {e.g., 
as parameterized by A^cff). However, some sterile neutrino parameters can lead to an acceptable 
A^eff ■ These parameter ranges are accompanied by considerable dilution of the ordinary background 
relic neutrinos, possibly an adverse effect on BBN, but sometimes fall in a range which can explain 
measured neutrino masses in some particle physics models. A robust signature of these sterile 
neutrinos would be a measured A^off 7^ 3 coupled with no cosmological signal for neutrino rest 
mass when the detection thresholds for these probes are below laboratory-established neutrino mass 
values, either as established by the atmospheric neutrino oscillation scale or direct measurements 
with, e.g., KATRIN or neutrino-less double beta decay experiments. 

PACS numbers; 14.60.Pq; 14.60.St; 26.35. +c; 95.30.-k; 95.30.Cq 



I. INTRODUCTION 

In this paper we show that recent Cosmic Microwave 
Background (CMB) limits on the proportion of energy 
density in the early universe contributed by particles with 
relativistic kinematics can provide powerful constraints 
on a possible sector of particle physics which is difh- 
cult and sometimes impossible to probe in the laboratory. 
The constraints we derive are independent of Big Bang 
Nucleosynthesis (BBN) considerations, though these may 
ultimately extend and strengthen our limits. 

Sterile neutrino effects on the conditions in BBN have 
been considered in Ref.s [HQ). Likewise, much work has 
been done on the effects of heavy particle decay on BBN 
and on particle decay-generated, post-BBN cascade nu- 
clear reactions and nucleosynthesis [sl-fT^. Here we con- 
centrate on a range of sterile neutrino masses and life- 
times which are, in some cases, different than what has 
been studied before and we concentrate on effects on the 
thermal neutrino background and the thermodynamics 
of the early universe. 

The experimental and observational establishment of 
neutrino flavor mixing and nonzero neutrino rest masses 
invites speculation about new beyond-Standard-Model 
physics in the neutrino sector. In particular, the exis- 
tence of right-handed neutral fermions (sterile neutrinos) 
seems at least plausible. Right-handed neutrinos with 
masses below the electroweak scale are called sterile neu- 
trinos. The LEP results require only three active neu- 



trinos with standard weak interactions, but there is no 
limit on the number of sterile neutrinos. There is no com- 
pelling argument for the rest mass scales of these states, 
with models invoking masses from the sub-eV range to 
the unification scale. Sterile neutrinos have been invoked 
as a way to engineer successful r-Process nucleosynthesis 
in neutrino- heated supernova ejecta [isl - ITTj . Sterile neu- 
trinos with keV rest masses have been studied as po- 
tential dark matter candidates I18ri33 , and as enablers of 
the formation of the first stars [3J-[3g , large pulsar kicks 
[37I - I41I , baryogenesis (i^ l43j , core c ollap se supernova ex- 
plosions (44| - |47t . and other issues [43]. (See a review 
on sterile neutrinos in Ref. (1^.) The results of some 
accelerator-based experiments have been interpreted as 
suggesting active-sterile neutrino flavor mixing [50l - l60j . 



The seesaw mechanism [61|-|65j can explain small active 
neutrino masses by postulating order one Yukawa cou- 
plings to the Higgs and very large right-handed neutrino 
masses. However, the split seesaw mechanism [66| can 
reconcile the small active neutrino masses and a sterile 
neutrino with a mass below the electroweak scale. Ideas 
on what might be the rest masses of sterile neutrinos 
abound [67,. ,68i] . 

For example, one class of models for generating the ob- 
served active neutrino masses (69| invokes heavy sterile 
neutrinos to induce a contribution to the lighter neutrino 
mass eigenvalues {i.e., those most closely associated with 
the active neutrinos). Interestingly, in our calculations of 
radiation energy density we find that the ranges of ster- 
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ile neutrino rest mass and vacuum mixing (with active 
species) which give the best agreement with observation 
sometimes fah in the "sweet spot" ranges which explain 
the observed neutrino masses in these models. 

Active-sterile neutrino mixing results in effective sub- 
weak interaction coupling between sterile neutrinos and 
the rest of the universe Jai\. The upshot is that the 
sterile states may be thermally populated in the early 
universe, and these sterile neutrinos can decay to leptons 
and photons. The key issue is that sterile neutrinos which 
mix with active species can play havoc with astrophysical 
environments where neutrinos influence energetics and 
composition. A case in point is the early universe. 

There are at least three classic observational avenues 
to probe or constrain a generation of particles in the 
early universe that have decayed away. First, the 
observationally-inferred light element abundances, when 
confronted with standard BBN calculations which uti- 
lize the CMB anisotropy-inferred baryon-to-photon ratio, 
could reveal discordance with observationally-inferred 
primordial light element abundances. Such discordance 
could, in turn, lead to constraint. 

Second, the decay products of these particles may not 
completely thermalize, leading to extra radiation energy 
density not accounted for by photons and the three ac- 
tive neutrino flavors. For example, WMAP7 reports 
iVoff = 4.34j;Q gg, where the radiation energy density prad 
is related to the photon temperature T through 
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A recent analysis of CMB and large scale structure data 
argues that iVoff = 4.08j:°:g8 El- Note that N^s is sim- 
ply a parametrization of the radiation energy density. In 
principle this parameter might have nothing to do with 
neutrinos, though the choice of factors in its definition 
is tied to the standard cosmology scenario in which it 
does. For example, the factor of {T^/Tf'^ = (4/11)^^/3 -^^ 
this definition comes from the standard cosmology case: 
Fermi-Dirac-shaped energy distribution functions for the 
ordinary relic neutrino backgrounds with "temperature" 
Ti, and degeneracy parameter zero. In the standard case 
A'cff should be close to 3, matching the number of active 
neutrino flavors. (A careful calculation in the standard 
cosmology case, accounting for de-coupling induced de- 
viations from thermal distribution functions for the neu- 
trinos, but retaining the definition of N^s in Eq.[l] yields 
iVeff = 3.046 ^.) 

Third, an epoch of matter domination by these parti- 
cles in the early universe could alter the growth or damp- 
ing of structures on scales smaller than that of the causal 
horizon at this epoch. Of course, this is of little use as 
a means of constraining the particles considered here be- 
cause the mass-scale of fluctuations so modified would 
be quite small, well below what can be probed via the 
CMB and current observations of small-scale, large-scale 
structure, e.g., the Lyman alpha forest. 



Here we assess constraints on sterile neutrino rest mass 
and lifetime parameter space mostly from the second of 
these considerations. To do this we sought to answer a 
seemingly straightforward question: What would be the 
effects on the thermal history of the early universe stem- 
ming from heavy sterile neutrinos with number densities 
of order those of photons but which decay at or around 
the BBN epoch in the early universe? We targeted ster- 
ile neutrinos with rest masses (100-500 Me V) very large 
compared to the BBN energy scale (~ IMeV), with 
vacuum mixing with active neutrinos sufficiently small 
that these species evade current laboratory neutrino mass 
bounds and drop out of thermal and chemical equilibrium 
at epochs with temperature T > 1 GeV, but sufficiently 
large that they would decay away during the BBN epoch. 
It is likely that the gross violence such a scenario would 
do to the time-temperature and scale factor-temperature 
relations and neutron-to-proton ratio during BBN would 
render the light element abundances sufficiently differ- 
ent from measured values so as to constrain this largely 
unexplored sector of neutrino physics. 

Even without light element BBN considerations, the 
work reported here demonstrates that a large range of the 
sterile neutrino rest mass and mixing parameter space, 
unconstrained by existing laboratory experiments, can 
nevertheless be ruled out via CMB limits on iVeff. 

Curiously, however, we find that there are ranges of 
these sterile neutrino parameters where NeS is consistent 
with observational constraints, but where the usual cos- 
mic relic active neutrino background makes a negligible 
contribution to A'cff and, instead, the dominant contri- 
bution to this quantity arises from the active neutrino 
decay products of the sterile neutrinos. We leave to a 
later paper whether scenarios with these sterile neutrino 
parameters can evade BBN-based primordial light ele- 
ment abundance bounds. 

Scenarios along these lines with decaying sterile neu- 
trinos that influence the relic neutrino background open 
a fourth channel for constraint: laboratory measurement 
of a neutrino rest mass and the comparison of this with 
the CMB- and structure-derived cosmological bounds on 
neutrino mass. A robust signature of the sterile neu- 
trino decay scenario outlined here would be a measured 
A'cff ^ 3 and a measurement of neutrino mass in, e.g., 
KATRIN or a neutrino-less double beta decay exper- 
iment, exceeding the cosmological bound on neutrino 
mass. 

In fact, it is better than this because we already know 
(from the square-root of the atmospheric neutrino mixing 
mass-squared splitting) that at least one neutrino mass 
eigenvalue exceeds 50 meV in the normal neutrino mass 
hierarchy, and two in the inverted hierarchy. As cos- 
mological neutrino mass constraints push down toward 
100 meV [tI] the stage is set for a definitive test of the 
entropy generation/dilution scenarios discussed here. 

We identify roughly the regions of sterile neutrino rest 
mass and decay lifetime (i.e., coupling or mixing with ac- 
tive neutrinos) that produce matter dominated regimes 
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that may last many Hubble times and where entropy 
production can be significant. On this latter point, the 
mechanism of out of equilibrium sterile neutrino decay- 
generated entropy production and its phasing with weak 
freeze-out and BBN is interesting in itself. 

In what follows we address current bounds on sterile 
neutrino masses and vacuum flavor mixing in section II, 
sterile neutrino decoupling in section III, entropy produc- 
tion and associated dilution from out-of-equilibrium ster- 
ile neutrino decay in section IV, iVcff modification result- 
ing from these decay scenarios and constraints in section 
V, speculation on the significance of "hidden" matter- 
dominated epochs in section VI, and dilution effects on 
the cosmological determination of active neutrino masses 
in section VII. We give conclusions in section VIII. 



II. CONSTRAINTS ON STERILE NEUTRINO 
REST MASSES AND FLAVOR MIXING 

A. Active neutrino mass limits 

Sterile neutrinos which mix in vacuum with active neu- 
trinos are not truly sterile because they have an effective 
interaction with strength ~ sin^ 9 Gp , where 9 is the 
effective vacuum mixing angle and the Fermi constant 
squared G| gives the strength of the ordinary weak in- 
teraction. This effective coupling, though perhaps very 
small, nonetheless provides an avenue for these particles 
to be created in the early universe, and constrained by 
direct laboratory experiments and observations. 

First consider the constraints. Laboratory neutrino 
mass constraints (tJ, [t^ , particularly those based on the 
tritium beta decay endpoint, which are good whether 
neutrinos are Majorana or Dirac in character, also con- 
strain the extent to which very massive sterile neutrinos 
can mix with, for example, the electron neutrino in vac- 
uum. The effective z/p "mass" in vacuum is 



[20|,l76|,l77 



\U, 



e4 



' 1714 = n^s sin^ 9, (2) 



where Ue i are the mass basis to flavor basis unitary trans- 
formation matrix elements, mi are the neutrino mass 
eigenvalues {i — 1,2,3,4) and we minimally assume a 
fourth of these, with the sterile neutrino mass roughly 
•nis = m4. Here we consider the limit where mi 3> toi^2.3, 
so that the effective vacuum active-sterile mixing an- 
gle is defined through sit? 9 = \Uei\^ ■ If m|j™ is the 
tritium endpoint experimental upper limit on neutrino 
mass, then we must have sin^ 9 < m]^"^/ms. 

Sterile neutrino masses less than the pion rest masses 
(m^o = 135 MeV, m^± = 139.57 MeV) wiU decay pri- 
marily through three neutrino decay Va + + Dp, 
while those above this mass-scale will decay via Vs — ?> 
T^'^ + Va, with a, /3 = e, /x, r, or via various charged pion 
decay modes that will be discussed below. The three- 
neutrino decay channel, for example, leads to lifetime 



4 ^MeVy 1 
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sm 



(3) 



while the tt*' decay channel, for example, produces an 
inverse decay rate which very roughly is 



r^o K. 5.76 X 10"^ s- 



1 



sm 



X {x^ — 1) ' 



(4) 



where x = mg/m^o. The other pion decay channels will 
be discussed below, but this one is adequate to illustrate 
rough laboratory TOs constraints. 

Neutrino mass upper limits from any laboratory source 
translate into lower bounds on the lifetimes of sterile neu- 
trinos, r: for the three neutrino decay mode. 



1 AA inio /'2eV\ /MeVV 
> 1.44 X 10^^ s ' ' ' ' ■ 



eV\ /MeVy 



while for the tt'^ decay mode 



> 0.389 s 



2eV\ 



Tlim 



- 1 



(5) 



(6) 



Neither of these constraints affects the sterile neutrino 
mass and lifetime parameter space of most interest in 
this paper. Accelerator-based laboratory constraints may 
impact the lowest mass and highest mixing angle sterile 
neutrinos considered here 'Ys'-SsJ. Future tritium end- 
point experiments like KATRIN may take the v^. mass 
limit down by an order of mag nitude, to m'™ w 0.2 eV, 
with appropriately more stringent lower limits on sterile 
neutrino lifetimes. Double beta decay experiments may 
provide far better limits, but these apply only to Majo- 
rana neutrinos. 

Cosmological constraints (based on large scale struc- 
ture and the CMB) on the sum of the light neutrino 
masses [tI, [s^ - lssj will come into our considerations in 
a manner very different from the usual arguments in 
standard cosmology. The active neutrinos decouple early 
(~ Is) and the usual cosmological constraints are fig- 
ured at photon decoupling {t ~ 10^^ s) or later, and at 
these epochs the sum of the light neutrino masses M^, 
will have little contribution from heavy sterile states, as 
these components will have decayed away. At time t we 
would have 



|C/,4|'e-*/" mi, 



(7) 



a—e,fi,T i— 1.2,3 



if the background neutrinos were in flavor states satis- 
fying all the usual cosmological assumptions and if mass 
state 4 mixes only with in vacuum, with obvious gener- 
alizations for more complicated active-sterile mixing ma- 
trix realizations. 

There are stringent cosmological constraints on long- 
lived extra sterile neutrino states, see e.g., Ref. [89[. 
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For the most part, these cosmological constraints are 
predicated on the assumptions of reUc active neutrino 
backgrounds with black body, thermal-shaped, Fermi- 
Dirac energy spectrum characterized by a temperature 
some 40% lower than the photon temperature, and zero 
chemical potentials. As we will show below, the sterile 
neutrino-decay scenarios considered here produce a very 
different relic active neutrino background. 



B. Neutrino Mass Models 

In some models for the origin of neutrino rest mass the 
ultimate mass eigenvalue is a sum of an intrinsic piece 
and a contribution induced by the presence of a mas- 
sive (mostly) sterile neutrino component as described in 
Ref. [6y| . This induced mass contribution is ^ nis 9^ , so 
that there are three broad categories we can consider. 

First is the regime where this product is much less 
than the known active neutrino mass scale, i.e., nis 0^ ^ 
0.05 eV, or where sin^ 261 < 10-^°/ (m^/lOO MeV). These 
neutrino mass models give effectively no constraint in this 
regime. Many if not most of the sterile neutrino mass and 
mixing/lifetime parameters we consider below fall in this 
regime, especially for nis > 150 MeV 

Second is the regime where 0.05 eV < ms9^ < 0.5 eV, 
or sin2 26l - (10-^" - 10"^) / (m^/lOO MeV). This is 
right in the range to give the requisite induced mass con- 
tribution to the active neutrinos. As will will see, this 
range actually overlaps with some of the ranges of ster- 
ile neutrino rest mass and mixing which give the best 
agreement with observational constraints on A^off- 

Finally, a third regime is where 9^ 3> 0.5 eV, corre- 
sponding to sin^ 29 > 10'^/ (mJlOOMeV). Sterile neu- 
trino rest mass and lifetime/mixing parameters falling in 
this range would be disfavored in these induced neutrino 
mass models. 



C. Considerations from Colliders 

Since the sterile neutrino has non-zero mixing with the 
active neutrino flavors, any interaction that produces an 
active neutrino can instead produce a sterile neutrino, 
with a cross section reduced by the effective vacuum 
active-sterile mixing, sin^ 9. If a sterile neutrino is cre- 
ated in a collider and then subsequently decays in a de- 
tector, a signature for the sterile neutrino may be in- 
ferred [905. For example, the decay of sterile neutrinos 
with mass larger than m^± will produce a co-linear pion- 
electron pair. 

The rate of sterile neutrino decays in the detector can 
be estimated from the overall sterile neutrino production 
rate by taking the product of the collider luminosity, L, 
the weak cross section, a ^ G^E^ , and the vacuum mix- 
ing angle suppression factor, sin^ 9. The detection rate 
would be the production rate, suppressed by the ratio of 
the crossing time of the neutrino through the detector 



to the laboratory-frame lifetime of the sterile neutrino 
(product of Lorentz factor 7 and t) 



R 

■JT 



R 



30 m / V 30 GeV 



V200 MeV/ VlOO sJ 



(8) 



Detection or constraints along these lines obtained from 
colliders would be difficult, but not impossible. 



III. STERILE NEUTRINO DECOUPLING IN 
THE EARLY UNIVERSE 

The sterile neutrinos discussed above can have effec- 
tive interactions large enough to keep them in thermal 
and chemical equilibrium for epochs in the early universe 
with temperatures sufficiently high. For the rest mass 
ranges we consider here these species will have relativistic 
kinematics when they decouple. Therefore, subsequently, 
they will have a momentum space distribution function 
that is a relativistic Fermi-Dirac black body, for which 
we will take the chemical potential to be zero. Of course, 
as the universe expands, coupled particles will annihilate 
or decrease in number and the entropy that these parti- 
cles carry will be transferred to the still-coupled particles 
and not to the decoupled sterile neutrinos. The net result 
will be that at the epoch of standard weak decoupling, 
temperature T ~ 1 MeV, the sterile neutrinos will have 
a "temperature" lower than that of the active neutrinos 
and the plasma. 

At early epochs, when the sterile neutrinos are coupled 
and possess relativistic kinematics, their scattering rate 
will be ~ G| sin^ 29. Thermal equilibrium between 
the sterile neutrinos and the plasma will obtain where 
this rate is larger than the local Hubble expansion rate. 
In radiation dominated conditions the expansion rate is 

H = (87r3/90)^^^ rVmpi, where g is the statisti- 
cal weight in relativistic particles and nipi is the Planck 
mass. Comparing these rates, the temperature where 
sterile neutrinos decouple is roughly 



^ dec 



87r3 
"90" 



2 GeV 



1/6 



9 



1/6 



V61.75/ 



(G| nipif^ sin2/3 29 
1/6 / in-9 \ 



sin^ 261 



(9) 
(10) 



For example, a sterile neutrino with rest mass rus — 
200 MeV and sin^ 26* = 10"^° would have lifetime t w 
43 s, and between its decoupling at T '-^ 4 GeV, and weak 
(active neutrino) decoupling at T w 3 MeV, a tenth of a 
second or so would have elapsed, and few sterile neu- 
trinos would have decayed. Instead, though they would 
then have largely nonrelativistic kinematics, their energy 
spectrum would have a relativistic Fermi-Dirac form, al- 
beit with a temperature lower than that for the plasma 
and active neutrinos. 
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Where particle decay is negligible the entropy in a co- 
moving volume will be conserved. If the entropy is carried 
by relativistic particles the ratio of the sterile neutrino 
temperature to the active neutrino and plasma tempera- 
tures at weak decoupling (wd) will be 



T 



wd 



1/3 



10.75 
61.75 



1/3 



1 



1.79' 



(11) 



where Qv^doc and ^wd are the statistical weight in rel- 
ativistic particles at the epochs of sterile neutrino and 
weak decoupling, respectively. The former quantity is 
figured to be roughly g^^dec « g-iefiv + 5qg: where the 
statistical weight in photons, neutrinos, e^, and is 
g-ye^v = 14.25, while that in quark and gluon degrees of 
freedom is gqg = 16-1- 10.5 iVf — 47.5, where we leave 
out tau leptons in the former and take three relativistic 
quark flavors, A^f = 3 {u,d,s), in the latter weight. By 
far the biggest source of dilution in T^_^ /T stems from the 
loss of quark and gluon degrees of freedom at the QCD 
epoch at T sa 170 MeV. 

Since the ordinary thermal background neutrinos have 
roughly the same temperature as photons at the begin- 
ning of weak decoupling. Eg. I 111 implies that the ratio of 
sterile neutrino to active background neutrino tempera- 
ture is 



G = 



T 



wd 



(12) 



Since the decoupled sterile and ordinary thermal back- 
ground neutrino energies redshift with scale factor in the 
same way, G will be a co-moving invariant, i.e., fixed. 

Despite dilution from epochs of particle annihilation, 
at weak decoupling the number density of sterile neutri- 
nos, assuming negligible numbers of decays, will still be 
a significant fraction of the photon number density n-y. 



3C(3)^3 



0.1 n, 



wd 



7' 



(13) 



where the Riemann Zeta function is ^(3) ~ 1.20206 and 
we include both left- and right-handed sterile neutrinos. 
As time goes on these sterile neutrinos will decay away. 
The actual proper number density of these particles at 
time t is 



3C(3) 



(14) 



Of course, there will be further dilution of T^JT as 
pairs annihilate and entropy is added to the pho- 
ton/baryon/electron plasma through sterile neutrino de- 
cay. 



IV. DECAY-INDUCED ENTROPY 
GENERATION AND DILUTION 

Heavy sterile neutrinos can decay out of equilibrium 
in the early universe and thereby generate entropy and 



cause what is sometimes termed dilution. These thermo- 
dynamic and cosmological issues were discussed in the 
general case for any unstable particles in Ref. [9l|. Here 
we examine dilution for the specific case of heavy sterile 
neutrinos decaying around the BBN epoch, but there are 
new twists in this case. For example, there are sterile 
neutrino decay channels in which ordinary active neutri- 
nos are produced. These decay-produced active neutri- 
nos may or may not thermalize. If they thermalize they 
add to dilution. If they do not they add to the sea of 
decoupled relativistic particles. Taken together, dilution 
and the addition of decoupled decay neutrinos affect A'off 
in competing ways. 



A. Sterile Neutrino Decay Processes 

Massive sterile neutrinos can decay through flavor mix- 
ing with active species as outlined above. If the sterile 
neutrino mass is less than the tt" mass, then these parti- 
cles will decay into three active neutrinos, with a smaller 
branch into a neutrino and a photon, or into a neutrino 
and an e^-pair. For sterile neutrinos with rest masses 
above the 7r° or rest mass scales other, more com- 
plicated decay channels beconie possible. These are dis- 
cussed below and in Ref.s and [tI,!!!]. 

For each decay mechanism, a fraction /em of the 
total sterile neutrino decay energy will be comprised 
of electromagnetic decay products (photons, electrons 
and positrons). This fraction of the decay energy will 
be added to the plasma through inelastic scattering of 
these products on particles in the photon-electron-baryon 
plasma. This is because the scattering timescale for these 
particles is much shorter than the dynamical timescale. 

On the other hand, the energy contained in the ac- 
tive neutrino decay products will remain in the neutrino 
sector, unless the individual neutrino energies are large 
enough that their scattering timescales rival the dynami- 
cal timescale. In the calculations presented here we com- 
pare the decay-produced neutrino scattering rate to the 
local Hubble expansion rate and estimate the neutrino 
energy where these rates are equal. 



rcthcrm 



-TT p 



1/2 



rp4 



(15) 



where p is the total mass-energy density, Gp is the Fermi 
constant, and mpi is the Planck mass. 

If a decay neutrino has an energy larger than £;rcthcrm 
we assume that it will down-scatter until it has energy 
^rothcrm_ rpj^^ difference between its initial energy and 
^rothorm .^jjj deposited in the plasma. The fraction of 
the sterile neutrino rest mass deposited in the plasma via 
the scattering of decay neutrinos is fi, , with / = /cm -I- f„ . 
Decay neutrinos with energies smaller than i<;™therm g^j.^ 
taken to be decoupled. Neutrinos which are decoupled 
do not contribute to the thermalization fraction /, but 
they do contribute to relativistic energy density, i.e., to 
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There are seven decay processes included in the calcu- 
lations presented here. These are: 



The decay rate for this process is 



192 7r3 



nCi ■ sin^ ( 



(16) 



3.47X lO-^s-i. f-^) -sm^e. (17) 



process, 



(21) 



The ratio of 1 /3 between the rates of two processes comes 
from the difference in the decay of the virtual in these 
cases: the decay is into three flavors of neutrino pairs in 
3-neutrino decay; and into a single e^-pair in the Vs — > 
1/ + e~ + e+ channel. We will have roughly f 2/3 for 
this process when rus is large enough. 



We make the assumption here that the heavy state 
decays through a coupling with one active species with 
effective vacuum mixing angle 9. This is tantamount to 
a particular mass basis/flavor basis transformation, with 
implications for degeneracy in the decay rate. The fi- 
nal state active decay neutrinos produced in this process 
usually are decoupled and so deposit no energy in the 
photon-electron-baryon plasma, i.e., / = 0. The excep- 
tion is for the most massive sterile neutrinos considered 
here, and only then very early in the expansion. 



2. 



1^ + "/ 



Here we follow Ref.s [20] and [77] and take the rate for 
this process as 



Gp c 

64 7r4 4 



_9G|_ 

512 TT* 



5 • 2 , 

a ■ TO„ • sm I 



2.72 X lO^'^s^i 



MeV 



sm 



(18) 

(19) 
(20) 



where a is the fine structure constant and we assume 
a particular vacuum flavor mixing structure, and again 
use rris = m^. This electromagnetic decay process is 
not GIM-suppressed because it involves sterile neutrinos 
with no corresponding charged lepton. In these expres- 
sions, F (r^) « — I -I- jrp, « {itlib/Mw)'^ , the ratio of 
charged lepton mass to the W mass, all squared, and the 
sum is over flavor /3. Because the sterile neutrinos con- 
sidered here have non-relativistic kinematics at relevant 
epochs, while the decay products are relativistic, the fi- 
nal state neutrino and photon each have energy ms/2. 
This implies that f Ki 1/2 for this process. 



Similarly, when the sterile neutrino rest mass is much 
larger than twice the muon rest mass (tos ^ 2 m^) the 
rate for this decay process is the same as for the last 
process. 



(22) 



Obviously this decay channel has a significant sterile neu- 
trino rest mass threshold, nis > 2m^ « 211.32 MeV. 
Moreover, the plasma energy deposition in this decay 
channel is different because the muons decay. For ex- 
ample, /i^ —5- + i/^ + i/e produces an electron and 
two neutrinos. The electron or positron thermalizes 
while the neutrinos, each of which has average energy 
Eu^^ = 34.33 MeV, usually remain decoupled and so do 
not. 



Sterile neutrinos with masses larger than the tt^ mass 
(tos > TO^o K, 135 MeV) can decay into a iv^ and an ac- 
tive neutrino. The tt" then decays into two photons. This 
decay channel has mean lifetime 8.4 x 10"^'' s, which is 
much shorter than both the scattering timescale and dy- 
namical timescale in the early universe, so we can regard 
these decays as being instantaneous. The energy of the 
active neutrino is 



ml - rnio 



2m., 



(23) 



The remaining available decay energy, m^. — E'^ , is ther- 
malized into the plasma as the tt'^ decays into photons. 
The decay rate for this channel is 



When the sterile neutrino rest mass is much larger than 
twice the electron rest mass, then the rate for this decay 
process is just 1/3 of the rate for the 3-neutrino decay 



r.o. - • [ml - mlo] ■ sin^ 9 (24) 

w 1.735 X lO^s^^a; [a;^ - 1] • sin^e*, (25) 

with X = ms/rriT^o and here we take /^r — 131 MeV. 
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FIG. 1: Branching ratios for selected decay processes and the 
branching ratio- weighted fraction of sterile neutrino rest mass 
(neglecting thermalization) deposited in the decoupled active 
neutrino seas are shown as functions of sterile neutrino rest 
mass in MeV. 



Sterile neutrinos with masses larger than the sum of 
the charged pion and electron masses (m^ > to^± + me ~ 
140.08 MeV) can decay into a tt^ and an e^. The electron 
(or positron) has energy 



2m, 



(26) 



which is quickly thermalized into the plasma. The 
charged pion then decays, e.g., 7r+ — ^ /i+ + (or its 
CP-conjugated counterpart), with a mean lifetime of 
2.6 X 10~*s, which is much shorter than the scattering 
and dynamical timescales, so we can also regard these de- 
cays as effectively instantaneous. The muon then decays 
producing two neutrinos as described above. Altogether 
three neutrinos are produced in this decay. Charged pion 
decay at rest produces a mono-energetic (or v^) with 
energy E^t^ = 29.79 MeV. The average energy added to 
the decoupled neutrino sea per decay in this channel is 
taken to be {Edcc} = E^.^ + 2Ey^ + (5/6)(£;^ke), where 
the kinetic energy of the pion is Et^ke = rUs — Eg— m^± . 



J 



The decay rate for this channel is 



16 TT 

1.917 X 10^ s 



2 / \2 



■(i£v)M('-<-'')00-<-»)0 



sm 



1/2 



sm 



(27) 
(28) 



with X + y = {rnT^± + mg) jnis and x — y = (m^i — me) /m^ and we again take = 131 MeV. Since there are two 
channels for this decay, i.e., into tt"*" or tt^, we take the total decay rate as F^e = 2F^+g-. 



The kinematics of this decay are similar to those of 
mass threshold, which is now m, > m^± -f m 



TT 



± 



-e^, but with the obvious difference in sterile neutrino 



p ~ 245.23MeV. The decay rate in this channel is 

1/2 



167r 

1.917 X 10® s"^ ( 



{rnl - (m^± -|- m^)^^ {rnl ~ (m,r± - m^)^ 
(l-{x + yf) (l~{x^y) 



MeV 



1/2 



sin^ e. 



(29) 
(30) 



with x + y = (jnT^± + m^) /m^ and x — y = (m^± — m^) /m^ and we again take /^r ~ 131 MeV. Since there are two 
channels for this decay, i.e., into 7r+ or tt", we take the total decay rate as F^^ — 2F^+„- . 



r 



The muon (or anti-muon) produced promptly in this 
decay has energy 



E„ 



2ms 



(31) 



The charged pion co-produced with this muon has to- 
tal energy Ej^ = m^ — Ef^ and decays producing another 
muon, e.g., 7r+ — > /i+ + (or its CP-conjugated coun- 
terpart). Subsequently the two muons decay, each pro- 
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ducing two neutrinos as described above. Altogether five 
neutrinos are produced in this decay sequence. The av- 
erage energy added to the decoupled neutrino sea per de- 
cay in this channel is taken to be (i?doc) = E^t^ + 4i?^^ + 
(5/6)(i?^KE) + (2/3)i?^KE, where the kinetic energy of 
the pion is Ej^ke = E-^ — 771^^+ , while the kinetic energy 
of the prompt muon is E^ke ~ Efj_ — m^. 

A caveat: the list of sterile neutrino decay processes 
presented here is not definitive nor is it exhaustive. For 
example, the sterile-active vacuum flavor mixing scheme 
remains a complete guess, and this guess affects degen- 
eracies in the rate calculations. (Do all active neutrino 
couple to the sterile in the same way, so that appropri- 
ate decay rates are multiplied by three?) Additionally, 
for sterile neutrino masses rus > 500 MeV, five-neutrino 
decay becomes dominant (2^ and even three-neutrino de- 
cay, with its five powers of nis dependence, will eventually 
overtake the pion decay channels with their three-powers 
dependence. Nevertheless, for the sterile rest mass ranges 
we consider here, where Nee and dilution effects are most 
interesting, the seven decay channels we employ suffice 
to rough out the key behavior. 

Figure [1] shows the branching ratios and the fraction 
of sterile neutrino energy deposited in the decoupled ac- 
tive neutrino sea all as functions of sterile neutrino rest 
mass. This figure does not show the (usually small) frac- 
tion of decay active neutrino energy which is thermalized. 
The sterile neutrino mass thresholds produce the obvious 
kings and ledges visible in this figure. For sterile neutri- 
nos with masses below the 7r° threshold, three-neutrino 
decay dominates and only roughly / « 17% of the rest 
mass ends up thermalized in the plasma, and this princi- 
pally through — > + + e+ and — > + 7. Above 
the 7r° and + thresholds it is a different story, 
with / ~ 60% until the sterile neutrino mass exceeds 
the TT^ -I- fT^ threshold where, because this decay channel 
produces five final state active neutrinos, / falls to less 
than 50%. Of course, this also means that the majority 
of energy deposited in the decoupled active neutrino sea 
comes from this process when rUs is high enough. 

Rates must be summed for all allowed decay channels 
to produce a total decay rate, i.e., an inverse lifetime l/r. 
The relation between the overall sterile neutrino decay 
lifetime r and the rest mass mg and sin^ 29 is shown in 
Fig. [5] and Fig. |31 (In these figures we plot sin^ 29 instead 
of sin^ 9 for ease of comparison to several extant sterile 
neutrino constraint works.) 



B. Entropy Generation 

The relation between the entropy-per-baryon S, in 
units of Boltzmann's constant fcb, and the baryon-to- 




. 2 
Sin 29 



FIG. 2: Contours of sterile neutrino lifetime r in seconds as 
functions of sterile neutrino rest mass nis and sin^ 29. 
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FIG. 3: Same as Fig. [H but for lower ranges of sterile neutrino 
rest mass and larger vacuum mixing angles. 
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photon ratio, rj = nt/nj, is 
^ 9s 



45C(3); 77 



(5.895 X 10' 



9\ (9. 



6.11 X 10" 



(32) 
(33) 



where Qs is the statistical weight in relativistic par- 
ticles carrying the entropy, so that S = S/ni, = 
(27r^/45)gsT'^/nb, where nb and are the proper num- 
ber densities of baryons and photons, respectively. In 
a standard universe at temperatures low enough that 
baryon number is conserved and with no particle de- 
cay, nuclear reactions, shocks, etc., s would be constant. 
Note, however, that even in this limit the relation be- 
tween 77 and s is not constant because gs changes with 
time/temperature. The Wagoner-Fowler-Hoyle h param- 
eter is defined as h = nb/{NpJ'^), where Np^ is Avo- 
gadro's number and Tg = T/IG^K. With this definition 
the entropy-per-baryon is s w 1.213 x lQ^{gs/2)/h and 
h = (3.368 X lO'' g cm"'^) ry, and the baryon rest mass 
density is pb[gcva~'^] = Tg h. 

Decay of decoupled, out of equilibrium sterile neutrinos 
will result in entropy being added to the plasma. The 
rate at which entropy-per-baryon is added is 



ds 
dt 



T 



,-(i-to)/T 



(34) 



where / is the fraction of the sterile neutrino decay en- 
ergy (the rest mass rus) which thermalizes in the plasma. 
In Eq. p4|) the number of sterile neutrinos {vg + i^s) per 
baryon at time to is xq = / A){T^JT)q/ rjQ, where the 
zero subscript means these quantities are evaluated at 
to, the age of the universe at which the calculation is 
started. In the calculations to follow we take Iq suffi- 
ciently early (e.g., well above the weak decoupling tem- 
perature) that prior decay can be neglected, in which 
case Xq sa (3/4)(l/1.79)3/ryo. 

Sterile neutrino decay adds entropy and so decreases 
the baryon-to-photon ratio, but in the end, at temper- 
atures far below those characteristic of the BBN epoch, 
we must get the CMB-determined values of these quan- 
tities, e. g., 7?wMAP ~ 6.11 X 10"^" as measured by 
WMAP (used to scale Eq. |33] above), implying 

swMAP ~ 5.895 X 10^. Therefore, scenarios in which 
entropy is added through particle decay must start out 
with a higher value of 770 (lower value of sq)- 

We have solved Eq. [Ml with a modified early universe 
expansion and Big Bang Nucleosynthesis code which 
treats all thermodynamic variables self consistently with 
the Friedman equation and all relevant weak interaction 
and sterile neutrino decay processes. 

Figure 2] shows the temperature of the plasma and the 
decoupled active neutrino seas as functions of decoupled 
active neutrino temperature, for several different scenar- 
ios, as calculated with our code. One scenario shown in 
this figure is just the standard radiation-dominated case. 
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FIG. 4: Temperature versus decoupled active neutrino tem- 
perature for several scenarios as labeled. The decoupled ac- 
tive neutrino temperature (dashed line), decreasing to the left 
here, is inversely proportional to scale factor, which therefore 
increases to the right. The dash-dot-dot (blue) line shows the 
standard expansion with no sterile neutrinos, exhibiting the 
transfer of entropy from electron/positron pairs to photons 
as the former annihilate. The lighter solid (red) curve shows 
what happens in a scenario with a sterile neutrino with rest 
mass nis = 275 MeV and lifetime r — 100 s. 



where annihilation occurs as electromagnetic equilib- 
rium shifts with expansion, transferring entropy to the 
photons, but not to the active neutrinos. This results in 
the apparent shallowing of the slope of the plasma (pho- 
ton) temperature near T ^ 100 keV. Entropy generated 
by sterile neutrino decay results in a similar, albeit more 
dramatic, phenomenon as is evident in Fig. SI where we 
show a scenario with particular sterile neutrino rest mass 
mc = 275 MeV and lifetime r = 100 s. 



C. Dilution 

It is evident in the sterile neutrino decay scenario de- 
picted in Fig. |4] that the fossil thermal neutrino relic 
background winds up considerably colder relative to the 
photons than in the standard cosmology case. This is 
dilution, a direct result of entropy generation, and anal- 
ogous to what happens in the post-inflation re-heating 
epoch. In the cases we consider here the entropy is car- 
ried by relativistic particles. The ratio of the after-to- 
before entropies-per-baryon we designate as the dilution 
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FIG. 5: Entropy-per-baryon (in units of Boltzmann's con- 
stant fcb) versus photon (plasma) temperature for the stan- 
dard cosmology case (constant co-moving entropy, dashed 
line) and for a scenario with a sterile neutrino with rest mass 
rua = 275 MeV and lifetime r = 100s (solid line). Begin- 
ning and ending entropy-per-baryon for the latter case as in- 
dicated. 
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FIG. 6; Contours (as labeled) of dilution factor F (the ratio 
of final-to-initial entropy-per-baryon) are given as functions 
of sterile neutrino rest mass in MeV and sin^ 26, where 6 is 
the characteristic effective two-by-two vacuum mixing angle 
between active neutrino species and the sterile species. 



factor 



F 



final 



•^initial 9si Ct? 



(35) 



where the indices i and f indicate initial (prior to any ster- 
ile neutrino decay) and final (long after sterile neutrino 
decay) values of statistical weight in relativistic particles 
gs, photon temperature T, and scale factor a, respec- 
tively. Standard cosmology, with a fixed co-moving en- 
tropy, will have F = \. The relationship between scale 
factor and temperature for decoupled particles {e.g., neu- 
trinos) is T^fttf — T^iUi, and we assume that the initial 
photon and neutrino temperatures are the same, so that 
the ratio of the final neutrino temperature to the final 
photon temperature is 



where F — 1, we recover the usual relationship be- 
tween the background neutrino temperature and the 
photon/plasma temperature at, e.g., photon decoupling, 

T^{/T{ = (4/11)* « 1/1.4, i.e., the final neutrino tem- 
perature should be 40% lower than the photon temper- 
ature. But if entropy has been generated by out-of- 
equilibrium particle decay we will have F > 1, with a con- 
sequently lower ratio of background neutrino-to-photon 
temperatures. Therefore, for example, the ratio of the 
neutrino temperature in the diluted scenario at the pho- 
ton decoupling epoch to the neutrino temperature at pho- 
ton decoupling in the standard cosmology {F = 1) case 
is 



dec 

^ V 

T17 dec Stan 
V 



(38) 



7k 

Tf 



1 

1 



9si 



1 



(36) 



(37) 



where in the last line we take as an example initial 
and final temperatures which bracket the BBN and e^- 
annihilation epochs. For the standard cosmology case. 



The photon decoupfing epoch is when the CMB photon 
temperature is T « 0.2 eV. 

Figure [5] shows the history of entropy addition for a 
particular case (sterile neutrino rest mass TOs = 275 MeV 
and lifetime r = 100 s) and for a standard constant co- 
moving entropy cosmology. For this case F = (5.9 x 
10'')/(3.3 X 10**) « 18, with F* w 2.6, which implies a 
ratio of final neutrino and photon temperatures 2.6 times 
lower than the standard case, roughly 1/3.6. 
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sin^ze 

FIG. 7: Contours (as labeled) of dilution factor F (the ratio 
of final-to-initial entropy-per-baryon) are given as functions 
of sterile neutrino rest mass in MeV and sin^ 26, where 9 is 
the characteristic effective two-by-two vacuum mixing angle 
between active neutrino species and the sterile species. 

In general longer sterile neutrino decay lifetime for a 
given rus results in more entropy generation, necessitat- 
ing lower starting entropy, hence higher starting 779. For 
some choices of sterile neutrino parameters, sterile decay 
can be responsible for 90% or more of swmap- Note that 
for the parameters chosen in the example in Fig. [SJ it is 
evident that the bulk of the entropy is added during the 
BBN epoch. This is true for a fair fraction of the sterile 
neutrino decay parameters considered here. BBN effects 
will be discussed in a subsequent paper. 



We can use the code discussed above to survey cos- 
mological effects for ranges of sterile neutrino masses 
and lifetimes. Dilution factors for a portions of these 
ranges are shown as contours in Fig IH] and Fig. [T] One 
trend in these figures is clear: the later the epoch at 
which the sterile neutrinos decay, the lower will be the 
temperature and, hence, the higher the ratio TOs/T and 
so the higher the entropy added per sterile neutrino de- 
cay. For sterile neutrinos with masses rrig < m^o three- 
neutrino decay will dominate the overall decay rate, im- 
plying relatively small / and, hence, modest dilution. 
Once the sterile neutrino rest masses are larger than the 
TOttO ~ 135 MeV mass threshold, but smaller than the 
efficient decay neutrino-generating processes which op- 
erate above the mT^± + nie ~ 140 MeV threshold, we see 
prodigious entropy generation and dilution. There is also 
substantial dilution and decay neutrino generation for 
longer lifetime sterile neutrinos with masses above the 
m^± + m./. « 245.23 MeV threshold. 

In Fig. [7] we see that the largest amounts of dilution 
correspond to ranges of sterile neutrino rest mass and 
vacuum mixing that fall in the "sweet spot" ranges which 
in some particle physics models [6^ explain the observed 
light neutrino mass scale. This is the second regime of 
TOs 0^ discussed in section II and, as we shall see below, 
this range also may give an acceptable A^off- 



V. GENERATION OF iVeft AND AN ALTERED 
RELIC NEUTRINO BACKGROUND 

Armed with the dilution factors and the average ac- 
tive neutrino energy per sterile neutrino decay, {E^°'^ (i)), 
injected into the decoupled seas of active neutrinos at 
time t, we can calculate the radiation energy density 
and, hence, iVes at the photon decoupling epoch. If 
at time t the thermal background neutrino and sterile 
neutrino temperatures are Ti, (t) and T^^ (i), respectively, 
and TJ is the thermal neutrino background tempera- 
ture at photon decoupling, then we can use the definition 
of iVoff in Eq. [T] to show that in time step At the incre- 
ment to iVeff is 



AiVeff(t) = ^ — ' ^ (39) 



120 l^^" 



dec Stan 



120 ( {Ef" it)) \ ^3 C (3) At \ ( nAt) \ ' . ( ^ 



77r2 V Tr^it) / \2 71-2 T J [t) J \^y7dccstan 



The first line in this equation is written in a way designed expression is a product of the neutrino energy injected 
to elucidate its physical meaning. The numerator in this per sterile neutrino decay and the second term, the num- 
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FIG. 8: Contours of Ncft as functions of sterile neutrino rest 
mass (in MeV), and sin'^ 26 as in Fig. [6] 



FIG. 9: Contours of Ncft as functions of sterile neutrino rest 
mass (in MeV), and sin^ 26 as in Fig. [8] but for a lower mass 
and larger vacuum mixing angle range. 



ber of decays per unit volume in At (number density of 
sterile species at time t multiplied by decay rate 1 /r and 
At), and a third term which is simply a redshift factor 
(energy density redshifts like four inverse powers of scale 
factor and scale factor is inversely proportional to decou- 
pled neutrino temperature). This numerator is therefore 
the energy density in active decay neutrinos produced by 
sterile neutrino decay at an earlier epoch t, redshifted 
down to the photon decoupling epoch. The factor in the 
denominator follows from the definition of NeS in Eq. [1] 
The second line, Eq. 1401 simply rearranges the terms 
in the first line to isolate co-moving invariants edec (t) = 
{Et^^{t))/T,{t) and G = (T,^ {t)lT,{t)). Using these 
co-moving quantities and Eq. [35] in Eq. HHl we can con- 
clude that 



eff 



dt 



180 C (3) 



— r • edcc (i) 
F i 



(41) 



We can use these considerations in the code described 
above to integrate this equation and derive N^s as a func- 
tion of sterile neutrino rest mass and lifetime. 

Fig. [8] and Fig.|9]show the results of these calculations. 
They give contours of iVoff as functions of sterile neutrino 
rest mass and sin^ 20, essentially in ranges of these pa- 
rameters identical to those used in the lifetime and dilu- 
tion factor figures above. These results are remarkable. 
Where dilution is large, yet not accompanied by signif- 
icant active decay neutrino generation {e.g., where mg 
is less than the charged pion decay channel thresholds. 



yet bigger than the tt" threshold), iVcff is smaller than 3, 
reflecting the dilution-caused reduction/refrigeration of 
the ordinary thermal neutrino background contribution 
to the radiation energy density at photon decoupling. On 
the other hand, for TOs large enough for sterile neutrino 
decay to proceed through — > tt^ -|- /x^, we will have 
prodigious active decay neutrino production and dilution- 
caused reduction of the the ordinary thermal background 
neutrino contribution to A'off. However, as is evident 
in Fig. [51 the former effect wins out over dilution and 
large values of A^cft can be generated, reflecting the large 
amount of decay neutrino energy injected into the active 
neutrino seas during the epoch of sterile neutrino decay 
in this case. 

The particular scenario with TOs — 275 MeV and r = 
100 s shown in Fig. 0] and Fig. [5] serves to illustrate these 
points. This case has dilution factor F — 17.6 and results 
in iVoff = 4.266, and of that only TV^'^crm ^ 0.0655 is con- 
tributed by the ordinary thermal background neutrinos. 
(This small contribution from the thermal background 
neutrinos follows from dilution: N^^™ w 3 • F^'^/'^ w 

3 • (17.6)""^^ « 0.0655.) 

It is a disturbing thought that there could be ranges 
of sterile neutrino rest mass and lifetime, unconstrained 
by laboratory experiments and existing cosmological 
bounds, with A^eff ~ 3 in seeming accord with standard 
cosmology, but where the ordinary thermal background 
neutrinos are grossly diluted, comprising only ^ 2% of 
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FIG. 10: The rehc neutrino energy spectrum for the particular case with = 275 MeV and r — 100 s. Here the neutrino 
number density per co-moving scaled neutrino energy (driv/de in MeV'^) is shown as a function of the dimensionless co-moving 
scaled neutrino energy parameter e — E^/T^ for three components: the ordinary thermal background neutrinos (short-dashed, 
blue line); sterile neutrino decay-generated neutrinos from the charged pion decay channels (solid, red line), and the tt" decay 
channel (dashed, green line). 



NcS- Perhaps BBN considerations will rule out these 
sterile neutrino parameters, but we already have con- 
straints on A^eff from CMB considerations (WMAP7 [TOj 
reports TVeff = 4.34^° j^) which promise to get much bet- 
ter with future observations. The current constraints, 
while not very good, already serve to eliminate much of 



Given the spectral energy distributions of the active 
neutrinos resulting from the various sterile neutrino de- 
cay channels [H, [olj, we can use the above considera- 
tions and calculate the actual spectrum of active neutri- 
nos left over from sterile neutrino decay. Fig. [TU] shows 
the the relic neutrino energy spectrum for the particu- 
lar case with rus = 275 MeV and t = 100 s. This fig- 
ure shows three components to the relic neutrino back- 
ground in this scenario: the ordinary thermal background 
{e.g., see Ref.s 95, 96] ); and the decay-generated neutri- 
nos from the charged pion decay channels and the tt*^ 
decay channel. For the sterile neutrino lifetime in this 
case, most of the decays will be where the age of the 
universe is order t > 100 s, where the average thermal 
neutrino energy will be < 100 keV, but the average 
decay neutrino energy will be more like 50 to 100 MeV 
in this case, roughly a thousand times the thermal neu- 
trino energy. By the time of photon decoupling, the ratio 



the sterile neutrino parameter space depicted in Fig. [8] 
and Fig. [HI especially for the higher and lower sterile 
neutrino rest masses on these figures and especially for 
relatively longer lifetimes (lower effective vacuum active- 
sterile mixing angles) on these plots. 



I 

of the thermal background neutrino temperature to the 
photon temperature in this case will be y7dcc^2n7doc _ 
(4/ll)(i/3).i^-i/3 ^ (4/11)1/3.(1/17.6)1/3 ^ Q Since 
rp-f dec ^ Q 2 qY this implies that the average energy of the 
thermal background neutrinos wiU be 3TJ'^°'^ - 0.16 eV, 
barely large enough to assure that these species are rela- 
tivistic at this epoch. The minimum neutrino mass eigen- 
value is the square root of the atmospheric neutrino mass- 
squared difference, y/Srr^^ ~ 0.055 eV. (The spectrum 
shown in Fig. [10] remains correct even if the kinematics 
of the neutrinos becomes non-relativistic, so long as 
is interpreted as the neutrino spacelike momentum mag- 
nitude.) Note that the peak of the decay neutrino relic 
spectrum is about a factor of 1000 larger than that of 
the background thermal neutrinos. This implies that the 
decay neutrinos have energies 160 eV at the photon 
decoupling epoch. For a maximum neutrino mass eigen- 
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value \J i5to^(.jjj, these neutrinos would still be relativistic 
at the current epoch, as the ratio of 160/. 055 « 2900, but 
the redshift at photon decoupling is only z-^dcc ~ 1100. 
Note also that in this scenario the number density of 
neutrinos in the thermal background is only about 6% 
(i.e., XjF K. 0.06) of the number density in a standard 
cosmology. Moreover, the total number density of decay 
neutrinos in the relic background is ^ 10% of that in the 
thermal neutrino relic background. 

The upshot of the number densities of the relic neutri- 
nos being low compared to the standard cosmology case 
and the decay neutrinos being at very much higher ener- 
gies than in this standard case is that it is unlikely that 
cosmological neutrino mass measurements will detect any 
neutrino mass in these sterile neutrino scenarios, even if 
A^cff is within observational bounds. 

In fact, a "smoking gun" signal for a dilution and par- 
ticle decay scenario along the lines of that discussed here 
would be the measurement of a radiation energy den- 
sity at photon decoupling with iVcff 7^ 3 and yet no 
detection of a neutrino mass even when the observa- 
tional sensitivity to the sum of the light neutrino masses 
is pushed below known, laboratory mass limits, either 
1\fbm^^ ^ 100 meV in the inverted neutrino mass hier- 
archy, \fim^^ ^ 50 meV in the normal mass hierarchy, 
or below a mass measured directly by tritium endpoint 
experiments (e.5., KATRIN) or neutrino- less double beta 
decay experiments. Prospects for bettering sensitivity to 
the sum of the neutrino masses in future observations 
and analyses are discussed in Ref. [73| . 



VI. STERILE NEUTRINO-ENGENDERED 
MATTER DOMINATED EPOCHS IN THE 
EARLY UNIVERSE 

For most of the ranges of sterile neutrino rest mass 
and decay lifetime discussed above there will be signifi- 
cant periods of matter domination in the early universe 
before, during, or even after the BBN epoch. This is not 
surprising because the sterile neutrinos are nearly as nu- 
merous as photons and they decay at temperature scales 
much lower than their rest masses. 

Fig. [TT] shows the ratio of the energy density in non- 
relativistic matter (sterile neutrino rest mass) as a func- 
tion of photon temperature for the particular case with 
sterile neutrino rest mass rris — 275 MeV and lifetime 
T = 100 s. In this scenario we see that matter can domi- 
nate over radiation by nearly an order of magnitude for 
what turns out to be many Hubble times at the sterile 
neutrino decay epoch. The peak of entropy generation 
in this scenario is, of course, during the peak epoch of 
sterile neutrino decay, which will tend to be somewhat 
later than peak matter domination. Moreover, the scale- 
factor/time relation is considerably altered here. 

In the particular case shown in Fig. [TT] the matter- 
dominated epoch occurs just prior to, and in the early 
stages of BBN. In a standard BBN scenario a scale- 
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FIG. 11: Ratio of matter-to-radiation energy density as a 
function of plasma/photon temperature (in MeV) for the in- 
dicated sterile neutrino rest mass, ttIs, and lifetime r. 

factor/time relation characteristic of matter-domination 
produces light element abundance yields which look noth- 
ing like the observationally-inferred primordial abun- 
dances, especially for ^He. In these sterile neutrino decay 
scenarios, however, entropy is being added to the system 
and BBN remains close to nuclear statistical equilibrium 
for the lighter elements, so it is not a foregone conclusion 
that BBN considerations can rule out all of the parame- 
ter space of sterile neutrino mass and lifetime not already 
ruled out by observationally-derived iVcff constraints. We 
will discuss detailed nucleosynthesis calculations for these 
scenarios in a subsequent work. 

Such drastic changes to the time-temperature-scale 
factor relationships that accompany sterile neutrino de- 
cay scenarios also result in altered causal and dynamical 
effects in the early universe. In broad brush, more energy 
density from sterile neutrino rest mass results in faster 
expansion rates and larger causal horizon lengths at a 
given temperature. But matter-domination adds an ad- 
ditional twist, driving the total mass-energy inside the 
horizon toward the Jeans mass. 

Fig. [T2] shows the total mass-energy inside the causal 
horizon as a function of photon temperature for the par- 
ticular scenario with sterile neutrino rest mass vris = 
275 MeV and lifetime r = 100 s. A causal horizon mass 
calculated using only the thermal radiation energy den- 
sity is shown to illustrate the effect of sterile neutrino 
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FIG. 12: Causal horizon mass (in solar masses) as a func- 
tion of plasma/photon temperature (in MeV) for a scenario 
(solid line) with the indicated sterile neutrino rest mass, rris, 
and lifetime r, and the horizon mass that would be inferred 
from the radiation energy density alone (dash-dotted line). 
The dashed line gives the Jeans mass for this sterile neutrino 
scenario. 



rest mass. The central take-awy message of this figure 
is that the mass scales affected by these sterile neutrino 
scenarios are small, of order the horizon scale Mh ~ 10^ 
to lO"^ M0, or smaller. These scales are marginally below 
what can currently be probed with CMB anisotropy and 
large scale structure arguments. 

This figure also shows the instantaneous Jeans mass as 
a function of temperature in the plasma. The Jeans mass 

M] — ^TT TOpj j / (6 y/p) , where Cs is the sound speed, is 

derived by setting the sound crossing time across a scale 
equal to the gravitational collapse (free fall) time on that 
scale. Generally, mass scales larger than the Jeans mass 
are subject to collapse, and fluctuations inside the hori- 
zon can grow in amplitude in matter-dominated condi- 
tions. In a radiation-dominated Friedman universe the 
causal horizon mass win be A/h (3/87r)^/^ (l/c3)Mj, 
i.e., less than the Jeans mass, and fluctuation growth is 
suppressed. 

In these matter-dominated sterile neutrino scenarios 
the sound speed is still set by radiation and remains close 
to light speed, Cg ~ l/V^ c. This means that, at best, 
the horizon mass and Jeans mass nearly coincide. This is 
evident in Fig. [T^ We conclude that collapse of horizon- 
sized regions is unlikely because there is no way for such 



structures to dump their entropy and cool. At best, the 
effect of the matter-dominated epoch engendered by ster- 
ile neutrinos might be to peak-up power on sub-horizon 
scales. These scales are relatively small. 



VII. IMPLICATIONS AND CONCLUSIONS 

We have assessed the impact of massive decaying ster- 
ile neutrinos on the thermal and dynamical history of the 
early universe and we have used our calculations to make 
predictions of fossil radiation energy density in these sce- 
narios which can be used to constrain them. In particu- 
lar, we have examined sterile neutrinos with rest masses 
in the 100 McV to 500 MeV range with vacuum cou- 
pling/mixing with ordinary active neutrinos large enough 
to engender thermal and chemical equilibrium for these 
sterile species at very early epochs in the history of the 
universe {i.e., T > IGeV), but small enough that they 
decay with lifetimes of seconds to minutes. 

We have shown that the out-of-equilibrium decay of 
these sterile species can lead to prodigious entropy pro- 
duction, as well as the generation of a significant popu- 
lation of decoupled active neutrinos. In these scenarios 
90% or more of the observed entropy-per-baryon in the 
universe can be generated by these decays. This neces- 
sarily means that pre-existing decoupled background par- 
ticles, e.g., the thermal background neutrinos, will suffer 
dilution, i.e., reduction in their number densities and en- 
ergies beyond that produced in a standard cosmological 
expansion. 

This can engineer a seemingly bizarre result. We find, 
for example, that ranges of sterile neutrino rest masses 
and lifetimes could produce a radiation energy density (as 
parameterized by Ncn) at the photon decoupling epoch 
within current observational bounds, but where the or- 
dinary thermal background neutrinos make only a negli- 
gible contribution to this. 

This dilution has a robust signature: a measured 
Ncff 7^ 3 (either smaller than 3 or significantly larger 
than this) coupled with no cosmological signal for neu- 
trino rest mass, even when the detection thresholds for 
the sum of the light neutrino masses in these probes 
are below laboratory-measured neutrino mass values, ei- 
ther as established by the atmospheric neutrino oscilla- 
tion scale (^ 50 meV in the normal neutrino mass hi- 
erarchy, ~ 100 meV in the inverted) or direct measure- 
ments with, e.g., tritium endpoint experiments like KA- 
TRIN or neutrino-less double beta decay experiments. 
There would be no neutrino mass detected in these sce- 
narios because cosmological probes are predicated on the 
assumption of black body, Fermi-Dirac-shaped spectral 
energy distributions (with temperature parameter some 
40% less than the photon background temperature) for 
the relic background neutrinos and this would not be the 
case with the significant dilution attending the sterile 
neutrino models considered here. 

Of course, many of the sterile neutrino rest mass and 
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lifetime parameters considered here will affect light el- 
ement abundance yields in Big Bang Nucleosynthesis 
(BBN). Likely this can result in additional constraints 
beyond the energy density constraints presented here. 
What, however, would be the implications of sterile neu- 
trino with masses and lifetimes which somehow managed 
to evade these constraints? 

The answer is straightforward if somewhat unsettling. 
We have shown that fluctuations on scales of ^ 10^ Mq 
could be affected, albeit seemingly minimally. More sig- 
nificantly, the massive dilution which is a by-product of 
sterile neutrino decay in these scenarios has two effects. 
First, since in these sterile neutrino scenarios there can be 
more than an order of magnitude increase in its baryon- 
to-entropy ratio, the pressure on existing baryogenesis 
models would be significantly increased. This is espe- 
cially true for electroweak baryogenesis schemes [97l4lOOl | . 

Second, dilution can hide many sins. Light decoupled 
particles could be "hidden" from cosmological detection 
and constraint by dilution. For example, light sterile neu- 
trinos, possibly those suggested by recent experiments 
|50l - l60j and some of those long popular with supernova 
heavy element nucleosynthesis modelers [isj . but seem- 
ingly ruled out by cosmological bounds, now might be 
compatible with these bounds. The dilution schemes dis- 
cussed here might then be an alternative to attempts to 
reconcile evidence for a light sterile with existing bounds 
[l0ll - [T03 |. 

There are caveats. Though dilution could reduce the 
number densities and energies of these particles enough 
to evade existing cosmological bounds, this is only strictly 
true if their abundances are set prior to the entropy 
generation/dilution epoch. If these light sterile species 
had large enough vacuum mixing angles then oscillations 
could regenerate appreciable populations of these par- 



ticles from the relic backgrounds left after dilution, and 
then the cosmological constraints might apply again if the 
relic decay neutrinos and thermal background neutrinos 
had a high enough number density post entropy genera- 
tion - which they may not have for enough dilution and 
low enough light sterile neutrino mass. Obviously BBN 
constraints on dilution scenarios can be another severe 
constraint, as it is not yet settled whether any of the 
dilution schemes discussed here can be reconciled with 
light element abundance bounds. 

Curiously, some of the ranges of sterile neutrino rest 
mass and vacuum mixing with active species which we 
have found to give acceptable values of N^s, and consid- 
erable dilution, also fall in the ranges of the parameters 
rec|uired by models to explain the light neutrino masses 

Of late we have learned much from the laboratory 
about the properties of neutrinos. However, the now 
established fact that neutrinos have rest masses opens 
a Pandora's Box of possibili ties for new physics in the 
neutrino sector [69l Il04l - ll06j |. These possibilities, while 
highlighting our ignorance, nevertheless may make for a 
future golden age of neutrin o physics when the next gen- 
eration cosmological probes |l07l4lld | become available. 
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